Long-range interaction of two metastable rare-gas atoms 
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We present semiempirical calculations of long-range van der Waals interactions for two interacting 
■ metastable rare-gas atoms Ne through Xe. Dispersion coefficients C'e are obtained for homonuclear 

molecular potentials asymptotically connecting to the ns(3/2)2 + ns(3/2)2 atomic states. The 
estimated uncertainty of the calculated Ce dispersion coefficients is 4%. 

& : 

Motivated by cold-collision studies of metastable rare-gas atoms and prospects of achieving Bose-Einstein 

Condensation in these systems ||, we present calculations of long-range dispersion (van der Waals) coefficients for 
two atoms interacting in the ns(3/2)2 atomic states (n = 3 for Ne, n = 4 for Ar, n = 5 for Kr, and n = 6 for Xe). The 
metastable states have long lifetimes, 43 sec for Xe ||, decaying to the ground ^So state by a weak magnetic-quadrupole 
transition. With such a long lifetime the metastable atom behaves as an effective ground state in experiments. 
O ,' Compared to alkali-metal systems, an attractive feature of the noble gas atoms is the availability of isotopes with 
zero nuclear spin. The lack of hyperfme structure leads to a substantial simplification of molecular potentials, though 
some complexity arises due to the nonvanishing total electron angular momentum (J=2) of the metastable state. 
The anisotropy leads to fifteen distinct long-range molecular states connecting to the ns(3/2)2 + ns(3/2)2 asymptotic 
. configuration. 

Our theoretical treatment of long-range interactions is similar to recent high-precision calculations of van der Waals 
coefficients for alkali-metal atoms |Q. By using many-body methods and accurate experimental matrix elements for 
the principal transitions, leading dispersion coefficients Cq were determined to an accuracy better than 1% for Na, 
K, and Rb, and of 1% for Cs and 1.5% for Fr. The semiempirical values of Cq coefficients for metastable noble-gas 
• atoms obtained here have an estimated uncertainty of 4%. The approach relies on the determination of dynamic 
polarizability functions. To construct the polarizabihties we combine experimental lifetime (S-15 and energy data of 



the excited states with accurate semiempirical dynamic polarizabihties of the ground states of noble-gas atoms [L6j. 
The theoretical lifetimes and branching ratios [17[ are adjusted to reproduce the measured static polarizabihties [18|, 
^ . which are known with a 2% uncertainty. We estimated the additional small contributions within the Dirac-Hartree- 
Fock framework. The resulting polarizabihties satisfy the Thomas-Reiche-Kuhn oscillator strength sum rule. 

The Racah notation for atomic levels is used. The particle-hole states are labeled as n£(K)j or n£'(K)j, where n 
and I are the principal and the orbital angular momentum quantum numbers of the valence electron and K = J c + £, 
where J c is the angular momentum of the core. The primed configurations converge to a Rydberg series limit with 
a hole in the (n — l)pi/ 2 state, and the unprimed to a hole in the (n — l)p3/ 2 state. The manifold of the lowest ns 
valence states has four fine-structure states ns'(l/2)o,i and ns(3/2)i.2, and the lowest np manifold consists of ten 
states. We investigate here the molecular potentials asymptotically connecting to the ns(3/2)2 atomic states. 

We calculate the long-range molecular potentials in the framework of Rayleigh-Schroedinger perturbation theory. 



The basis functions are defined as products of atomic wavefunctions 

\M 1 M 2 ;fl) = |ns(3/2) 2 M 1 ) 1 |ns(3/2) 2 M 2 ) 2 , (1) 

where the index 1(2) describes the wavefunction located on the center 1(2) and fl = Mi + M 2 , Mi l2 being projections 
of the atomic total angular momentum on the internuclear axis. Due to the axial symmetry of a dimer 17 is a conserved 
quantum number, ft takes values ranging from zero to four. The two-atom basis ([!]) is degenerate and the correct 
molecular wavefunctions are obtained by diagonalizing the molecular Hamiltonian 

H = H 1 +H 2 + V(R) . (2) 

In expression (|^) Hk represent the Hamiltonians of the two non-interacting atoms, and V(R) is the interaction 
potential at an internuclear distance R. The energy of the ns(3/2)2 metastable state is designated as £ * . Then in the 
model space (fy) 



Hx + H 2 j \MiM 2 ; fi) = 2£* \MxM 2 ] . 

The residual electrostatic potential V(R) is defined as the full Coulomb interaction energy in the dimer excluding 
interactions of the atomic electrons with their parent nuclei. 
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The multipole interactions (L — 1 for dipole-dipole, and L = 2 for quadrupole-quadrupole interactions) are given 

by @ 



VMR) = ^ £ ( £ -gfl' + „), (*>) , • < 3 » 



with the multipole spherical tensors 



Tf^-N^rfcW^), (4) 



where the summation is over atomic electrons, is the position vector of electron i, and C^ L \ii) are reduced spherical 
harmonics [2C| ]. In the following we write = and Q M = t/ 2 -* 



The lowest-order contribution to the term energies arises from the quadrupole-quadrupole interaction V qqi which 
varies as 1/R 5 . However, the corresponding C5 coefficients are only of the order 10 1 — 10 2 a.u. ||, and the dominant 
contribution appears in the second order in V(R), arising from the dipole-dipole interaction Vdd- The second-order 
dipole interaction is proportional to 1/R 6 , and the associated dispersion coefficient Cq is of the order of 10 4 — 10 5 a.u.. 
Applying the formalism of degenerate perturbation theory in second order we obtain an effective Hamiltonian 
within the two-atom basis Eq. (|l|) 

H2$ In) = 2rS mn + (m\V qq \n) + £ ^f^f^ . (5) 

The intermediate molecular state with unperturbed energy Ei runs over a complete set of two-atom states, 
excluding the model-space states Eq. ([!]). The formalism of the generalized Bloch equation 0] would allow the 
inclusion in the model space of the other three atomic states in the ns manifold and would account for the mixing 
of the different fine-structure levels; but such a large model space is not necessary for R > 10 a.u.. The position of 
the avoided level crossing can be estimated from i?j. c . « 2Cg/ {£ n s(3/2) 1 — £ns(3/2) 2 ) ~ 10 a - u -- The fl = 4 molecular 
term is unique and being unaffected by avoided crossings, the region of applicability is extended to R ~ n a.u., before 
the electronic clouds start to overlap. The effect of the quadrupole-quadrupole interaction on the term energy can 
be disregarded at values of R <C Cq/Cs ~ 10 3 a.u.. The quadrupole-quadrupole correction is discussed by Doery et 
al @. 

Using the Wigner-Eckart theorem, we can represent the matrix element of the dipole-dipole term in the effective 
second-order Hamiltonian as 



E 



(MiM 2 ; »|T/ drf |^)(^|U dd |Af{A^; fl) = 
2£* - Ei 

-4 E c^i-iy^W ^H A t( M ^ M i) A -\-^ Q M ^ M 'i) • ( 6 ) 

II' X/jl 

The dipole weights u>* are w^_ 1 = w 1 l 1 = 1, and Wq = 2. J a and Jt are the corresponding total angular momenta of 
intermediate atomic states of atoms 1 and 2, and 



(7) 



where m = 5 (Mi + M[ + A — fi) . The intermediate (uncoupled) dispersion coefficients are 



c J.J b = _ duSj (iw)S Jb (iu). (8) 

7T JO 

The reduced dynamic dipole polarizability Sj{iu>) of purely imaginary argument is defined as the sum over atomic 
states \aiIMj) with total angular momentum / and energy £ ai 



Sl{lUj) = E (f*-f ai )=« +w 2 ■ ( 9 ) 
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Here aj stands for all quantum numbers of the intermediate state, except for the total angular momentum /, and 
are the reduced electric-dipole matrix elements, defined by the Wigner-Eckart theorem. Three sums with 
I = 1,2,3 are allowed by electric-dipole selection rules. 

We proceed now to construct the dipole polarizability functions Si(iu>), Eq. (O), evaluate the uncoupled dispersion 
coefficients (||), and set and diagonalize the second-order effective Hamiltonian (|5|) . 

The functions Si{iuj) satisfy several sum rules. In particular, the static tensor dipole polarizability a zz (M) of the 
ns(3/2)2 state may be expressed as 

a zz {M) = -2 £(-!)'( _m I * 5/(0). (10) 

The static tensor dipole polarizabilities a zz (M = 1) and a zz (M — 2) of metastable noble-gas atoms Ne through Xe 
have been measured by Molof et al. |18) to within an error of 2%. In the present calculations the values of 5/(0) 
are adjusted to reproduce these experimental values. In addition, as to — > oo, the reduced polarizabilities satisfy the 
nonrelativistic Thomas-Reiche-Kuhn (TRK) sum rule 

= ^J2(-l) I+1 Sdioo) =N, (11) 

n I 

N being the number of electrons in the atom. Our constructed polarizabilities satisfy the sum rule. 

It is instructive to consider the action of a one-particle operator on the reference particle-hole Slater determinant 
ns (3/2)2 in the independent electron approximation. Such an operator can (i) annihilate the reference particle- hole 
pair; (ii) promote a valence electron from ns\/2 state to another valence state TOPi/2,3/2) the state of the (n — l)p3/2 
hole being unchanged; (iii) deexcite the (n — l)f>3/2 hole into some other hole state, the valence state remaining the 
same, and (iv) create another particle-hole pair in addition to the reference pair. According to such a classification it 
is convenient to break the polarizability function, Eq. (0), into three contributions corresponding to a number k 
of particle-hole pairs in the intermediate state \ot-iI) 

5/ = S°j+ Sj + Sj . 

Since an electric-dipole transition from ns(3/2)2 to the closed-core state ^So is prohibited by the angular selection 
rules, S°j = 0. 

The sum Sj is separated into contributions from the intermediate states in the lowest np fine-structure multiplct 
(Sj) , and the rest of the sum (S})' 

s} = (s}) np + (s})'. 

The first term is the dominant contribution. We calculate [Sj) using experimental values of transition energies and 

decay rates, and adjusted branching ratios. The rest of the sum over valence states (including bound and continuum 
states) (Sj)' is estimated in the Dirac-Hartree-Fock (DHF) approximation. The metastable state \vh\J) in lowest 
order is represented as a combination of the \v) = nsi/ 2 particle state and the hole state \h) = (n — 1)^3/2, coupled to 
the total angular momentum J — 2, [((n — l)p3/2) _1 ^s 1 / 2 ]2- The lowest-order energy of such a state is £* = e v — £h, 
Ei being the energy of the DHF orbital The intermediate state is represented as \ma; J), a particle state m coupled 
with hole state a to the total momentum /. Explicitly, 



Pin*,) = (-D J+J [/][j]E{^ ( T £m) ?!"" )2 1 



+6. 



(e v -e m ) (v\\d\\m) 2 f J 1 I 

Jm Jh Jv 

(e a -e h )(/i||d||a> 2 f J 1 I ^ 2 



mv 



e h f +lu 2 1 3a iv jh 



where [K] = 2K + 1, the summation is performed over the core orbitals a and excited states m, excluding states of 
the lowest np multiplet, and J = 2. The first sum is associated with excitation of the valence electron (case (ii)) while 
the second sum with deexcitation of the hole state (case (iii)). To arrive at this result we disregarded the coupling 
between levels within the same fine-structure multiplet. For example, for Ne the [(2p3/ 2 ) _1 4pi/2]27 [(2p3/2) _1 4p3/2]2j 
and [(2pi/2)~ 1 4p3/2]2 states are summed over independently, even though the correct lowest-order wave-function is 
a linear combination of them. This approximation corresponds to a disregard of the small difference between the 
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energies of the coupled and uncoupled states. Since the contribution (Sj)' is relatively small, such an estimate suffices 
at the present level of accuracy. Numerical evaluation of (S})' has been performed using a B-spline basis set in the 
Vjsr-i DHF potential, with the hole in the (n — 1)^3/2 core orbital p3|. 

We separate the sum Sj over the core excited states into two contributions 



2\ 



sf = (sf) rnr +(sf 



cntr 




The first term is associated with the dynamic polarizability a g (iuj) of the closed-shell ground state ^So and the second 
term is a corrective counter term. 

(5 / 2 ) corc M = (-l) /+J+1 Ha s M. (12) 

We use the semiempirical dynamic polarizabilities for the ground states of noble-gas atoms of Kumar and Meath |L6| . 
The estimated uncertainty of these core polarizabilities is less than 1%. The primary role of the core polarizability 
is to provide the correct limit Eq. ( fj"l| ) at w — > 00. The relative importance of the core-excitation contribution 
increases for heavier systems; for example, in a similar calculation for Fr Q|, core excitations contribute 23% of the 
Cg dispersion coefficient. The high-frequency limit, Eq. ([ll]), is accurately reproduced by the present total reduced 
dynamic polarizabilities Si(iuj). We obtain for Ne 9.98, for Ar 17.95, for Kr 35.95, and for Xe 53.96 compared to the 
nonrelativistically exact values 10, 18, 36, and 54 respectively. 

In the ns(3/2)2 = [((n — l)P3/2)~ ln Si /2]2 state core excitations to the occupied magnetic sub-states of the ns\/2 
particle state are not allowed by the Pauli exclusion principle, and neither are the core excitations from the empty 
magnetic sub-state of the hole (n — 1)^3/2 ■ To remove these transitions from the core polarizability contribution 
(^)corc' we introduce a counter term (Sj) cntr - Explicitly in the independent-electron model 

(sf) t M = [i] [j] (-i) I+J ( V li J 1 

(fm - £h)(h\\d\\m) 2 

We estimate the small counter-term using the Dirac-Hartree-Fock approximation. 

The largest contribution to the sums Sj arises from the intermediate states in the lowest np fine-structure multiplet. 
The determination of electric-dipole matrix elements involved in the sum ($}) np requires a knowledge of both decay 
rates and branching ratios in the manifold. The relevant lifetimes have been measured to within an error less than 1% 
for Ne |]||, Ar [|l0| and Kr (ll]], and less than 3% for Xe |T^] . However, the branching ratios B are not established 
to the same precision. The most accurate measurements of B in Ne have an error bar of approximately 4-5%, 

which would introduce an uncertainty of 4-5% in the static polarizabilities, and 8-10% inaccuracy in the values of Cq. 
To reduce the consequent errors, the experimental values of the static polarizability, accurate to 2%, were chosen as 
the reference data. 

The branching ratios of transitions to the ns(3/2)2 state have been adjusted as follows. The sum 5*3(0) includes 
only one intermediate state in the np manifold, np(b/2)z, and very small ab initio corrections. The np{^>/2) 3 state 
has a single decay channel, so that the sum £3(0) is known with the experimental precision of the decay rate. The 
sums S 2 (0) and Si (0) can be deduced from the experimental values of the static tensor polarizability as 

5i(0) = ~S 3 (0) + 5a„(l) - I a zz (2) 
S 2 (0) = ^S 3 (0)-^a zz (2). 

Removing small ab initio and semiempirical core-excitation contributions from these sums, the sums (Si) (0) and 
(S2) np (0) are obtained. The branching ratios B for four states involved in the J = 1 sum and three states in 
the J — 2 sum were multiplied by a uniform scaling factor. Branching ratios for Ne || for the J = 1 levels were 
multiplied by 1.0035, and for the J = 2 level by 0.905 in order to reproduce the experimental values of the static 
tensor polarizabilities. We modified the recommended values of B for Ar ]15| ] by multiplying the branching ratios of 
the J = 1 states by 0.973 and of the J = 2 states by 0.965; the values used in the calculations are listed in Table 1. 
For Kr the velocity- gauge branching ratios, tabulated in Ref. from calculations by Aymar and Coulombe JL7| , 
were multiplied by 1.127 for the J = 1 states and by 1.0016 for the J = 2 states. For Xe, velocity-gauge values of B 
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calculated in Ref. (l7) were multiplied by 0.927 for the J = 1 states and by 0.929 for the J — 2 states. The adjusted 
data for Ar and Xe are listed in Table §. Doery et al. || have compiled the input data for Ne and Kr, which have to 
be similarly modified. 

We employ the constructed reduced polarizabilities Si(iu>) to calculate the intermediate uncoupled dispersion co- 
efficients CQ aJb by quadrature using Eq. (||). The coefficients are listed in Table ||. They are to be used if the entire 
molecular Hamiltonian, including quadrupole-quadrupole and higher multipoles or perturbation-theory orders is to 

be diagonalized. Finally, the molecular terms are obtained by the diagonalization of H^S , given by Eq. (g). Neglect 
of the small corrections due to the quadrupole-quadrupole interaction results in parameterization of term energies in 
the form 



U(R) = 2£* - Cq/R 6 



The calculated dispersion coefficients Cq for various molecular symmetries are listed in Table [II. Since the region 
close to uj = contributes the most to the values of the integral in Eq. (||), the uncertainty in the values of Cg is 
approximately 4%, reflecting the 2% experimental error in the static dipole tensor polarizabilities fl8| . The values 
of the Cq coefficients grow monotonically from Ne to Xe, due to the reduction in the energy separations between 
the metastable states and the np-manifold. For heavier systems the anisotropy in Cq, arising from relativistic effects 
becomes increasingly marked, from 6.5% in Ne to 16% in Xe. 

Long-range dispersion coefficients for two interacting metastable Ne atoms were evaluated recently by Doery et 
al. H . The Cq coefficients were calculated from the diagonalization of the molecular dipole-dipole Hamiltonian in the 
model space containing the lowest np-manifold, so limiting the intermediate states to the lowest np-manifold states 
in the present formulation. Experimental values of decay rates and branching ratios were used to deduce the electric- 
dipole matrix elements. The precision of the calculated values of Cq is about 8-10% due to the large uncertainty in 
the branching ratios. The values of Cq from Ref. || for different molecular symmetries vary between 1951 and 1956 
a.u., exhibiting much less anisotropy than the present results which range between 1877 and 1999 a.u.. The difference 
can be traced to the anisotropy in the static dipole polarizabilities. Indeed, utilizing input data from Ref. || we 
obtain a zz (M = 1) = 192 and a zz (M = 2) = 189 a.u. if we include only the np manifold as in Ref. ||. While 
ct zz (M = 1) agrees with the experimental value jL8| 192(4), the a zz (M = 2) is overestimated by three standard 
deviations compared to the experimental value 180(3) a.u.. 

The accuracy of the dispersion coefficients could be improved by applying relativistic all-order many-body meth- 
ods [^3 24 1 to calculate transition amplitudes between ns — np manifolds. Such ab initio calculations are intrinsically 
more challenging than for alkali- metal atoms; the accurate experimental lifetimes would provide an excellent gauge 
of accuracy. 

Our values of Cq coefficients will be useful in studies of cold collisions of metastable rare-gas atoms For 
example, we can estimate the rate coefficient for Penning ionization by ignoring spin-polarization and assuming that 
every trajectory that surmounts the angular momentum barrier leads to ionization p],p5|. The corresponding rate 
coefficient is given by p5fl 



k = 6.35 x 10~ 9 6 l/o cm 3 s" 1 . 

1/2 



where fi is the reduced mass measured in units of the electron mass and T is the temperature. Combined with 
short-range potentials pp^] a number of other properties could be determined. For example, scattering lengths of 
elastic collisions could be found, providing input for mean-field equations describing dilute quantum gases. 

This work was supported by the U.S. Department of Energy, Division of Chemical Sciences, Office of Energy 
Research. Thanks are due to M. R. Doery, S. Kotochigova, and J. F. Babb for useful discussions. The authors are 
grateful to W. R. Johnson for providing the B-spline routine for the V/v-i Dirac-Hartree-Fock potential. 



TABLE I. Input data for Ar and Xe calculations. Lifetimes r for Ar are from Volz and Schmoranzer, Ref. |10| , and for Xe 
from Inoue et al., Ref. (l2j, except where noted. Branching ratios B to the ns(3/2)2 level are adjusted to reproduce experimental 
tensor dipole polarizabilities, as discussed in the text. 



state 


r, ns 


Ar, n = 4 


B{%) 


Xe, n — 6 

t, ns 


B(%) 


np'(l/2)i 


27.85(7) 




17.66 


43.5(1. 5) b 


5.15 


np'(3/2) 2 


29.01(7) 




10.59 


38.1(1.3) c 


2.30 


np'(3/2)i 


29.83(8) 




1.84 


49(2) 


1.16 


np(3/2) 2 


28.52(7) 




68.55 


31(1) 


65.42 


np(3/2)i 


29.62(7) 




14.97 


37(1) 


7.45 
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np(5/2) 2 31.17(7) 27.76 39(1) 33.29 

np(5/2) 3 29.00(7) 100.00 31(1) 100.00 

np(l/2)i 39.2(2.2) a 72.11 38(1) 87.60 



a Wiese et at, Ref. |Q. 
b Allen et al, Ref. jfj. 
c Husson and Margerie, Ref. Jl3| ]. 



TABLE II. Intermediate dispersion coefficients Cg 7 " 76 , a.u., multiplied by a factor 10 



s-ill /^i2l s~i22 <^>31 /— *32 z^33 
^8 ^6 ^6 ^6 <-^6 

Ne 4.945 -7.333 10.89 10.92 -16.21 24.14 

Ar 12.48 -17.42 24.39 27.12 -37.88 59.02 

Kr 15.12 -20.18 27.01 31.80 -42.44 66.96 

Xe 23.24 -29.22 36.91 47.74 -60.06 98.31 



TABLE III. 


Dispersion coefficients Ce in a.u. 


for the interaction of two metastable 


1 ns(3/2)2 noble-gas atoms. 




Term 


Ne 


Ar 


Kr 


Xe 




1877 


4417 


4994 


7138 


3 S 


1919 


4565 


5195 


7490 


3<u 


1922 


4583 


5224 


7557 


2 9 


1967 


4751 


5459 


7991 


2 9 


1935 


4629 


5286 


7664 




1934 


4623 


5276 


7641 


i« 


1983 


4811 


5543 


8148 


i 9 


1982 


4810 


5541 


8145 


i 9 


1920 


4574 


5210 


7524 




1920 


4574 


5210 


7526 


°i 


1999 


4872 


5629 


8311 




1968 


4756 


5467 


8010 




1966 


4747 


5452 


7975 




1877 


4418 


4996 


7140 


o 9 + 


1877 


4418 


4996 


7140 
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